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Abstract 

Cosmic strings provide a radically different paradigm for the formation of structure 
to the prevailing inflationary one. They afford some extra technical complications: for 
example, the calculation of the power spectrum of matter and radiation perturbations 
requires the knowledge of the history of the evolution of the defects in the form of two- 
time correlation functions. We describe some numerical simulations of string networks, 
designed to measure the two-time correlations during their evolution. 



1 Introduction 



Calculations of the Cosmic Microwave Background (CMB) fluctuations from cosmic strings 
[m have proved much harder than their inflationary equivalents. One of the reasons is that 
the energy-momentum of the defects cannot be considered to be a perfect fluid like most of 
the rest of the contents of the Universe: instead we have to use numerical simulations to find 
how important quantities such as the energy density and velocity evolve. The defect energy- 
momentum is then a source for perturbations in the gravitational field, which must be calculated 
either using Greens functions or by direct numerical integration [^]. 

A major problem with using numerical simulations is the lack of dynamic range, even with 
supercomputers. The largest simulations can only span a ratio of conformal times of about 
10, which is not totally satisfactory. For calculating simple quantities such as power spectra 
an alternative approach is available: to use numerical simulations to construct an accurate 
model of the appropriate source functions, which are two-time correlation functions of various 
components of the energy-momentum tensor. Relatively simple Greens functions can then be 



used to compute the power spectra of perturbations in perfect fluids, including CDM and tightly 
coupled photons and baryons. 

This talk summarises some work in progress on the properties of cosmic strings as sources 
of the gravitational perturbations. An issue of interest which we address is that of coherence, 
recently raised in 0, and explored in more detail by the same authors in 0. 



2 Perturbations from defects 

A simple model serves to illustrate the issues involved. We suppose that the Universe is spa- 
tially fiat with zero cosmological constant, and consists of Cold Dark Matter (CDM), photons, 
baryons, and defects, with average densities pc, p-y, Pb and ps respectively. The fluid density 
perturbations are written 6c, 6^, 6b and the velocity perturbations are v^, v^, vb- The defect 
energy-momentum tensor is 0^,^, which is also a small perturbation to the background, and to 
first order can be considered separately conserved (or We take the synchronous gauge. 

Before decoupling at a redshift of z^ec — 1100 (assuming the standard ionization history), the 
photons and baryons were tightly coupled by Thompson scattering, which forces 6b = j6^ and 
Vb = v^. 

It is very convenient to use an entropy s = (|5^ — 6c), and its time derivative s, which is 
equal to the divergence of the radiation peculiar velocity. It is also convenient to introduce the 
pseudoenergy Too [0], which is part of an ordinarily conserved energy-momentum pseudotensor: 

where R = 3pB/4p^. In the sychronous gauge the pseudo-energy is in fact proportional to the 
scalar curvature of the constant conformal time slices. 
Let us now define the pseudo-energy perturbation 

We may use the conservation equations for the defects to write the equations of motion of the 
fluids as four first-order equations 

L = ^Q)[5oo-(^^c + 2fi^(l + i?))5,-2(]^(l + i?)s]-^(^^^f7A, (3) 
s = y, (4) 
= (^^{1 + ?>w)6qq + ^Vl^{1 + R)y + n,kv„ (5) 

ij = .k^cl{s + 6c)~{-]{l-^cl)y, (6) 



a 



where w = p/ p, the ratio of the total pressure to the total energy density, and = 1/3(1 + R). 

These equations have the inhomogeneous form Aj = M-Aj + Ej, where A = {6c, s, 6oo, y)^, 
and the source vector, picked out in bold above, is S = {—^{a/a)QsSs, 0, flgkvs, 0)^. The 
solution to these equations with initial condition A{r]i) is then 

A(r^) = G{v,v,)A{7^.) + r rfVG(r^, V)S(V), (7) 



where G{ri, rj') is the Green's function. The initial perturbation is A(?7j) = {6c{rii), 0, 0, 0)-^, cor- 
responding to perturbations which are both adiabatic, in the sense that s = 0, and isocurvature. 



in the sense that Too = 0. In the presence of defects, the isocurvature condition forces the initial 
perturbations in the matter and radiation to compensate those of the defects. However, the 
initial compensation can be ignored, as the initial condition consists only of decaying modes, 
which is one of the advantages of using this basis for the perturbations An initial condition 
consisting purely of adiabatic growing modes would be A(?7j) = ((5c(?7i), 0, 2(5c(?7j), 0)^. 
At time rj, the power spectra of the various perturbations are therefore 

= r rfr^i r dv,Gi{v,Vi)Gl{v,V2){^AVi)K{V2)), (8) 

Jrji Jrii 

(no sum on i). It is clear from this equation that the calculation of the power spectra requires 
the knowledge of the two-time correlation functions 

C^'{Vi,V2) = {6s{vimV2)), C^^(r]i,r/2) = {vs{viK{V2)), X''^(r]i, r^s) = {SsiviVM), 

(9) 

where the angle brackets are to be understood as averages over ensembles of topological defects. 
3 Flat space cosmic string simulations 

In order to compute the averages in (|^) we performed many simulations of networks of string 
in Minkowski space 0, |[. It is not at first sight such a good idea to abandon the expanding 
Universe, but there are arguments which support the procedure. Firstly, a spatially flat FRW 
cosmology is conformal to Minkowksi space, so one can identify Minkowksi time with conformal 
time, and Minkowksi space coordinates with comoving space coordinates. Secondly, the string 
network scale length is significantly smaller than the horizon length, which means that the 
effects of the expanding background on the dynamics of the string should be small. Thirdly, 
a string network can be well described by a couple of parameters: the string density and 
the correlation length of the network 0, which means we can hope to translate Minkowksi 
space results into FRW results by appropriate scalings. Lastly, the savings in computational 
resources are immense. The strings can be evolved on a cubic spatial lattice, which means that 
only integer arithmetic need be used, and that self-intersections are extremely easy to check 
for. 

A vital property possessed by string networks is known as scaling. In its simplest form, 
scaling means that all quantities with dimensions of length are proportional to a fundamental 
network scale ^, which is in turn a function of length scales in the network dynamics. In a 
cosmological setting, the only scale in the equations of motion is the horizon or the conformal 
time 1]. Thus c<. rj, where ^ is understood to be a comoving scale. The invariant length density 
L (which is proportional to the energy density through the string mass per unit length /i) is 
therefore, on naive dimensional grounds, proportional to In fact, it is convenient to define 
^ by ^ = l/^/Ti. The simulations conserve energy and thus invariant length. However, real 
strings decay into particles and gravitational radiation, so we model this by excluding small 
loops below a threshold Ic from being counted as string. This threshold may be chosen in 
many ways: it can be a constant, or it may be a fixed fraction of ^. However, we find that the 
remainder ("long" or "infinite" string) always scales towards the end of the simulations, that 
is, = XT], although the constant x depends on how the cut-off is chosen. 

Scaling has important implications for the correlation functions. For example, the correct 
scaling behaviour for the power spectrum of the string density is 



(10) 



Figure 1: Two-time correlation function (Boo (/c, ''7)600 ''?')) V — 225 



We can also express the velocity power spectrum in terms of a similar scaling function P^{kri). 
These forms ensure that the mean square fluctuation in the string density at horizon crossing 
is constant, when ^/rj is constant. Hence the compensating fluctuations in the matter and 
radiation are also constant at horizon crossing, which gives a Harrison-Zel'dovich spectrum. 

In Figures 1, 2, and 3 the correlation functions (|^) measured for strings whose length is 
above the threshold Ic are displayed. They result from many simulations on (645)^ or (1285)'^ 
lattices, where 6 is the lattice spacing, with approximately 7500 or 60000 string points. For 
calculating the two-time correlation functions we typically averaged over 50 simulations.^ We 
find that good fits in the important regions where the correlators are large are obtained from 
the following functions 



1 



r^k^{l-{kA))(r,-r,')'^ 



= ^pu{kOP^{k^')e- 



(11) 

(12) 
(13) 



The values for T, T' and T" are given in Table |l|. A is approximately 36, which indicates that 
it is probably a lattice artifact. 

Table 1. The values of the parameters for the models in equations (PTD, (|T^ and ([T^), 

obtained by minimising the x^- 



Ic 


T 


a 






25 


0.21 ±0.05 


0.19 ±0.05 


0.42 ± 0.05 


0.36 ±0.07 


45 


0.18 ±0.05 


0.16 ±0.05 


0.42 ± 0.05 


0.42 ±0.06 



4 Implications for strings in FRW backgrounds 

The effect of the expansion of the Universe in a Friedmann model is to reduce the velocity string 
segments by Hubble damping, just as for particles. However, as the string moves relativistically, 

^For more details of the simulations, the reader is referred to B. 





Figure 3: Two-time correlation function {U {k , ri)U* {k , rj')) for t]' = 226 

the equation of motion becomes non-linear. The conservation equations for Gqo and U = ihiQ^i 
become 



eoo + - (Ooo + e) + A;?7 = 0, 
a 

f7 + 2- f/ + -A;(e + 26^) = 0, 
a 3 



(14) 
(15) 



where 9 = Qu and 9** = {kikj — ^6ij)Qij. Both 9 and 9*^ are unconstrained by energy- 
momentum conservation, and so one could imagine that fluctuations in the pressure could 
drive extra fluctuations in the energy density through the second term in (|T^. As the scale of 
this term is r]~^ rather than k, this could change the conclusions of the previous section, that 
the coherence time for a mode of wavenumber k was proportional to k~^. Instead we could 
have the situation assumed in 0, with coherence time proportional to rj. Outside the horizon, 
this makes a significant difference. 

Let us examine how this might work, for superhorizon modes [ki] <^ 1). Firstly, we divide 
9oo and 9 into coherent and incoherent parts: 



9 



00 



9oo + A 



00, 



9 = SwsQoo + r. 



(16) 



where r is a random variable, whose fluctuations are limited only by the requirement of scale 
invar iance. Thus 



{r{rii)r*{ri2)) 



(17) 



It is then not hard to show that energy-momentum conservation imphes that 

(lAooP) - (|rp). 

Since (|rp) < (|0p), the extra fluctuations in the density induced by the pressure term are 
controlled by the size of the pressure term. Now, in the Minkowski space simulations, we 
find that (|ep) ~ 0.1(|eooP) for krj < 1, and thus we argue that any extra fluctuations in 
Friedmann models, with a coherence time set by the horizon rj, are likely to be small. 



5 Conclusions 



We have taken the first steps on the road to calculating the matter and radiation power spectra 
in the cosmic string scenario, by constructing a realistic model of the important parts of the 
string energy-momentum tensor which drive the fluid perturbations. We have argued that the 
Minkowski space simulations we used to construct the model incorporates the essential features 
of string evolution in an expanding universe, and that corrections due to the expansion are 
small. This of course should be checked with FRW string codes. Perhaps the most important 
feature of the string network is its coherence time, the time over which the phases of the 
Fourier components of the string energy-momentum tensor are correlated. We find that there 
is no single coherence time for the whole network. Although the fall-off away from the equal- 
time value is modulated by functions whose dominant behaviour is exp(— (r^i — ri2)'^/2ril), and 
the scale is set by for each mode, rjc is different for different correlation functions. 
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